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ABSTRACT 

We present direct spectroscopic modeling of 11 high-S/N observed spectra of the 
Type la SN 2010kg, taken between -10 and -1-5 days with respect to B-maximum. The 
synthetic spectra, calculated with the SYN-I—I- code, span the range between 4100 
and 8500 A. Our results are in good agreement with previous findings for other Type 
la SNe. Most of the spectral features are formed at or close to the photosphere, but 
some ions, like Fell and Mg 11 , also form features at ^^2000 - 5000 km s“^ above the 
photosphere. The well-known high-velocity features of the Call IR-triplet as well as 
Sill A6355 are also detected. 

The single absorption feature at ~4400 A, which usually has been identified as due 
to Si 111, is poorly fit with Si ill in SN 2010kg. We find that the fit can be improved by 
assuming that this feature is due to either C ill or O ii, located in the outermost part 
of the ejecta, ^4000 - 5000 km s“^ above the photosphere. Since the presence of Cm 
is unlikely, because of the lack of the necessary excitation/ionization conditions in the 
outer ejecta, we identify this feature as due to On. The simultaneous presence of Oi 
and O II is in good agreement with the optical depth calculations and the temperature 
distribution in the ejecta of SN 2010kg. This could be the first identification of singly 
ionized oxygen in a Type la SN atmosphere. 

Key words: (stars:) supernovae: individual: SN 2010kg - stars: abundances - line: 
identification - line: profiles 


1 INTRODUCTION 

The model of Type la supernovae (SNe la) is widely ac¬ 
cepted to be a thermonuclear explosion of a C-0 white 
dwarf (Hoyle & Fowler 1960), although there are several 
open questions about the progenitor system and the details 
of explosion. The cosmological application of la SNe is based 
on the standardization of their luminosity (Matheson et al. 
2012) because their peak luminosities are well correlated 
with the decline rates (Amis)!?)) of their light curves 
(Phillips et al. 1993). In order to increase the effectiveness 
of la SNe as distance indicators, we have to hnd the connec¬ 
tion between the observed diversity of la SNe (Branch et al. 
2006; Blondin et al. 2012) and the different explosion mod¬ 
els and/or combustion methods. The two frequently con¬ 
sidered explosion channels are the single-degenerate (the 
WD accretes material from its non-degenerate companion; 
Whelan & Iben 1973) and the double-degenerate (the pro¬ 
genitor merges with another WD; Iben & Tutukov 1984; 
Webbink 1984) model. The burning propagation can oc¬ 


cur as deflagration (propagation with subsonic speed; 
Nomoto et al. 1984) or delayed detonation (subsonic propa¬ 
gation turns into supersonic; Khokhlov 1991). The study of 
the early spectra and the investigation of the ion signatures 
may lead us to the explanation of the la SNe diversity. 


In this paper, we present the direct analysis of eleven 
spectra of the Type la SN 2010kg. After the overview of our 
dataset in Section 2, our spectrum htting method is pre¬ 
sented, focusing on the ion identification and the limitation 
of the approximations used by SYN-I—I- (Thomas et al. 2011) 
in Section 3. The results of the fitting process are presented 
and discussed in Section 4, while in Section 5 we focus on the 
temperature profile of SN 2010kg and the possible presence 
of unburnt material in the outer ejecta. Finally, we summa¬ 
rize our conclusions in Section 6. 
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Table 1. Journal of the spectroscopic observations. The columns contain the followings: date of observation, modified Julian date, phase 
with respect to the moment of maximum light in B-band, exposure time, spectral range, spectral resolution in A and signal-to-noise 
ratio, respectively. 


Date 

MJD 

Phase (days) 

Exp. time (s) 

Spectral range (A) 

AA (A) 

S/N 

2010-12-02 

55532.201 

-10 

1200 

4100-10,500 

19 

77 

2010-12-04 

55534.196 

-8 

1100 

4100-10,500 

19 

110 

2010-12-06 

55536.190 

-6 

900 

4100-10,500 

19 

108 

2010-12-07 

55537.190 

-5 

900 

4100-10,500 

19 

117 

2010-12-08 

55538.183 

-4 

700 

4100-10,500 

19 

98 

2010-12-09 

55539.193 

-3 

700 

4100-10,500 

19 

104 

2010-12-11 

55541.174 

-1 

2235 

4100-10,500 

19 

106 

2010-12-12 

55542.174 

0 

600 

4100-10,500 

19 

113 

2010-12-13 

55543.170 

+1 

600 

4100-10,500 

19 

79 

2010-12-15 

55545.164 

+3 

600 

4100-10,500 

19 

83 

2010-12-17 

55547.167 

+5 

600 

4100-10,500 

19 

90 


2 DATASET 3 METHODS 


SN 2010kg was discovered on 30th November 2010 
(Nayak et al. 2010), by the Lick Observatory Supernova 
Search team with the Katzman Automatic Imaging Tele¬ 
scope. The supernova exploded in NGC 1633, with redshift 
of z = 0.0166 (Springob et al. 2005). From the spectroscopic 
measurements, it was identified as a broad-line Type la su¬ 
pernova at more than 1 week before maximum brightness 
(Marion et al. 2010; Silverman et al. 2010). 

Broad-line (BL) SNe are a subclass of Type la, de¬ 
fined by Branch et al. (2006), based on the strength of the 
Sill A6355 A feature, which is much broader than in core¬ 
normal SNe la. As former studies showed (Blondin et al. 
2013, 2015), broad Sill A6355 lines in Type la SNe could 
be the result of delayed detonations. The spectra usually 
display high mean expansion velocities for all line-forming 
regions, while the decline rate of the velocities show large 
variation (Blondin et al. 2012; Silverman et al. 2012). 


Our data sample consists of eleven spectra of SN 2010kg 
between -10 and -|-5 days with respect to the moment of 
maximum luminosity in B-band (Dec 12.5, 2010). All of 
them were obtained with the Marcario Low Resolution Spec¬ 
trograph attached to the Hobby-Eberly Telescope (HET) at 
McDonald Observatory, Texas (Hill et al. 1998). All spec¬ 
tra extend from 4,100 to 10,500 A with resolution of ~19 
A and have signal-to-noise ratio between 70 and 120 (see 
Table 1). Before the analysis, all spectra were corrected for 
redshift and the Milky Way interstellar reddening, assuming 
Rv = 3.1 and the extinction curve of Fitzpatrick & Massa 
(2007), based on dust emission maps from COBE/DIRBE 
and IRAS/ISSA (Schlafly & Finkbeiner 2011). We applied 
E{B — V) = 0.137 mag for dereddening. Due to the cut¬ 
off at 4100 A in the HET spectra, we are able to analyse 
some of the strongest features, such as the Ca ll near-infrared 
triplet and Si ll A6355, but not the Ca ll H&K lines. Since the 
flux calibration becomes unreliable near the edges of the ob¬ 
served spectra, we restrict our analysis to the spectral range 
between 4300 and 8500 A. 


3.1 SYN-I-+ 

We use the SYN-I—I- code (Thomas et al. 2011), which is the 
modern version of the original SYNOW spectrum synthesis 
code (Fisher et al. 1997) rewritten in C-|—1-. The computa¬ 
tion follows the Sobolev-approximation (Sobolev 1960): the 
blackbody radiation emitted by the sharp, spherical photo¬ 
sphere interacts with the expanding atmosphere, which has 
a large velocity gradient and low density. The main sim¬ 
plification is that a photon interacts with an atom only 
in a narrow region within the ejecta, where it suffers res¬ 
onance scattering. The result is a P Cygni profile, where 
the pseudo-absorption is blueshifted with respect to the line 
central wavelength. 

The global input parameters for the whole spectrum 
model are the velocity (hereafter Vphot) and the tempera¬ 
ture (Tphot) of the photosphere. SYN-I—I- approximates the 
photosphere as a sharp, spherically symmetric surface, which 
emits pure blackbody flux and its blackbody temperature is 
Tphot- Although the shape of the SN spectra can be more-or- 
less described as a diluted blackbody (e.g. Jeffery & Branch 
1990), LTE conditions probably never occur in Type la 
ejecta, thus, the derived Tphot values may not represent the 
actual local temperatures (see also Sec. 5.1). 

Since SYN-I—I- computes relative fluxes, we have to mul¬ 
tiply the absolute fluxes of the observed spectra with a flux 
scaling parameter (aO in the input file of the code), which 
also needs to be set for each spectrum. SYN-I—I- is also able 
to scale the spectrum with a linear and quadratic function 
of the wavelength, but we set these warping parameters (al 
and a2 in the input file) to zero in order to reproduce the 
shape of the blackbody continuum without additional dis¬ 
tortion. There are additional parameters for every ion: the 
optical depth (r), the minimum velocity of the line-forming 
region (hereafter minimum velocity or Vmin), the e-folding 
velocity of the opacity profile (aux in the input file) and 
the Boltzmann excitation temperature (Texc)- We fix the 
maximum velocity for all line-forming regions at Vmax = 
40,000 km s“^. All these ion-parameters refer to a refer¬ 
ence line of the given ion. The optical depths of every other 
lines are calculated based on the T^xc parameter and the 
atomic data included by SYN-f-f. The applied atomic data 
were downloaded from the homepage of ES: Elementary Su- 
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pernova Spectrum Synthesis User Guide^, which is the en¬ 
hanced version of database of Kurucz (1993). 

The fitting code, SYNAPPS, varies the above men¬ 
tioned parameters automatically and simultaneously, and 
fits the model to the data via x^-minimisation. A caveat 
in such a complex fitting is that the large number of param¬ 
eters may make the final best-fit model ambiguous, because 
the fitting algorithm could stick in an incorrect parameter 
combination due to the strong correlations between the pa¬ 
rameters. With this in our view, we used mainly the SYN++ 
code instead of SYNAPPS and varied the input parameters 
interactively to get better fits. 


3.2 Ion identification 

During the photospheric phase, the P Cygni line profiles are 
dominant in the Type la SN spectrum. We focus on the spec¬ 
tra taken earlier than +14 days post-maximum, because the 
physical assumption for the line forming mechanism (pure 
resonant scattering) in the SYN++ code may be more-or- 
less valid roughly before this epoch. At the moment of maxi¬ 
mum light, the Call H&k lines at A3934 and A3969, the ’W” 
feature caused by Sll at ~5500 A, the Sill A6355 and the 
Call NIR triplet are always observable. 

Most of the lines in the supernova spectrum overlap 
each other due to the fast expansion of the ejecta, making the 
line identification always ambiguous. Owing to the several 
independent studies (e.g. Branch et al. 2006; Tanaka et al. 
2008; Parrent et al. 2012; Hsiao et al. 2013; Marion et al. 
2015; Jack et al. 2015), we can assume the existence of a 
few types of ions in the ejecta. These are typically Sill, 
Sll, Call (as mentioned above) supplemented by Fell, Felll, 
Mg II and O l. Some studies pointed out the possible contri¬ 
bution of Nal (Branch et al. 2005), Cl (Hsiao et al. 2015), 
Cll (Silverman & Filippenko 2012), Till (Filippenko et al. 
1992), Nill and Coll (Branch et al. 2005) at least in a frac¬ 
tion of la spectra taken before maximum light. 

We performed the ion identification based on the lit¬ 
erature mentioned above, and the detailed optical depth 
calculations presented by Hatano et al. (1999a). The reli¬ 
ably identified ions are Call, Sill, Sll, Fell, Felll, Ol and 
Mg II. Furthermore, we use Nal and Till for some synthetic 
spectra to get better fitting, but their existence - without a 
strong absorption peak - remains ambiguous. Contribution 
from Nill and Coll are also tested in order to reproduce 
the strong flux decrease below 4000 A observed in all Type 
la SNe (Branch & Venketakrishna 1986; Foley et al. 2012a). 
The contributions from all the individual ions at two differ¬ 
ent epochs are presented in Fig. 1 and 2. 

An open question concerning the line identification in 
Type la SNe is the origin of the deep, narrow absorption 
feature at ~4400 A. Looking through the literature, many 
authors fit this region with SI ill (e.g. Marion et al. 2014; 
Yamanaka et al. 2009) and with a mixed contribution of 
Coll and Fell. We find that, at least in the case of SN 
2010kg, either Cm or On may give a better fit than Sim; 
in Section 5.3 we discuss the possible absence / presence of 
Cm/ O II, respectively. 


^ https://c3.lbl.gov/es/ 


3.3 High Velocity Features 

High Velocity Features (HVF) have been the focus of interest 
since their first detection (Hatano et al. 1999b). Those ions, 
with Vmin being ^7000 - 15,000 km s“^ higher than the pho¬ 
tospheric velocity (vphot), produce more highly blueshifted 
absorption features (e.g. Silverman et al. 2015, hereafter 
S15; Mazzali et al. 2005). Together with the second compo¬ 
nent formed close to the photosphere (Photospheric Velocity 
Feature, PVF), they produce a typical double-bottom line 
profile. In homologously expanding ejecta, the different ve¬ 
locity values indicate that the line-forming regions of HVF 
and PVF are physically distinct, which is supported by the 
spectropolarimetric observations (Maund et al. 2013). The 
origin of HVFs are still uncertain; density or abundance en¬ 
hancement, caused by swept-up gas (Gerardy et al. 2004) or 
clumpy circumstellar material (CSM) (Kasen et al. 2003), 
could explain their origin. 

After fifteen years of research it seems certain that Call 
HVFs appear in most of the SNe la ejecta at early epochs 
and they become weaker as the SN evolves. The only ex¬ 
ception is the group of the peculiar 91bg-like SNe la where 
almost no HVFs are seen (S15). Sin HVFs are more rare. 
Approaching maximum light. Sill HVF usually disappear, 
and only Call HVF survives the maximum (Childress et al. 
2014). Some studies have also shown that other ions, like 
Fell (Marion et al. 2014) or Ol (Parrent et al. 2012), can 
produce HVFs a few days after the explosion. Unfortunately, 
the broad, overlapping lines can easily blend with these fea¬ 
tures, thus the identification of these HVFs in the case of 
other ions is still ambiguous in SN 2010kg. In our study, 
we investigate all the ions at higher velocities to create con¬ 
straints on the real HVFs. 

Nevertheless, the situation (and thus the whole spec¬ 
tral modelling) is more complex because of the possible 
presence of detached features (DF, Jeffery & Branch 1990). 
Their line-forming regions are above the photosphere, so the 
produced DFs show higher velocities than Vphot, but these 
velocities stay below the Vmin of HVFs. Thus, unlike the 
HVFs, the origin of DFs is unlikely to be related to the 
presence of CSM. The minimum velocity of a DF is at least 
~1500 - 2000 km s“^ higher than Vphot at any epochs and 
roughly lower than ~19 - 20,000 km s“^ one week before 
maximum light. Note that the upper velocity limit of the 
DF region depends on the epoch and cannot be fixed to 
an exact velocity value. Thus, we identify the features with 
Vphot < Vmin < Vphot + 2000 km s“^ as photospheric veloc¬ 
ity features (PVFs), while those having Vmin > Vphot + 6000 
km s“^ as HVFs. In between these two regimes, the features 
are considered as DFs. 

Since the origin of HVFs is uncertain, and they might be 
formed outside the SN ejecta (e.g. due to CSM-interaction), 
we refer the line forming region of the DFs with highest Vmin 
(with ~5 - 6000 km s“^ above the photosphere) as the outer 
region of the SN ejecta hereafter. Based on the recent SN 
la models (e.g. Dessart et al. 2014c; Seitenzahl et al. 2013), 
oxygen and carbon DFs may be expected, because of their 
higher mass fraction in the outer layers [vmin > 16,000 km 
s“^), , but these are also somewhat model dependent be¬ 
cause of the strength of mixing involved in the models (see 
e.g. Fig.2 of Dessart et al. 2014b). 
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o 

Rest wavelength (A) 


Figure 1. The single-ion synthetic spectra ten days before maximum light in SN 2010kg; the observed spectrum (grey line) and the 
whole synthetic spectrum (red line) 
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o 

Rest wavelength (A) 


Figure 2. The single-ion synthetic spectra at the maximum light in SN 2010kg; the observed spectrum (grey line) and the whole synthetic 
spectrum (red line) 
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Figure 3. The velocities of the high-velocity features and their photospheric components in SN 2010kg, compared to the results for SN 
2010kg of S15. The photospheric velocity is presented with grey solid line. 


4 RESULTS 

Adopting the velocity values from the fits, we are able 
roughly to map the spatial distribution of the ions in the 
ejecta. This can help us getting abundance information and 
creating constraints on the explosion mechanism. In the fol¬ 
lowing, we discuss each individual ion identified in the spec¬ 
tra of SN 2010kg. 


4.1 Ca II lines 

The singly-ionized calcium produces the strongest HVFs and 
the highest minimum velocities in Type la SNe spectra. Sev¬ 
eral studies (e.g. Childress et al. 2014; S15) showed that the 
Call HVFs can survive even after maximum light. These 
characteristics make the Ca ll H&K lines and the NIR triplet 
ideal for studying HVFs. The HET spectra used in this study 
do not cover the region below 4100 A, thus, Call H&K lines 
cannot be examined. The NIR triplet is surpassingly strong 
compared to spectra of other Type la SNe. The triplet is 
built up by Call A8498, A8542 and A8662. Although, be¬ 
cause of the high optical depth values needed at the first 
epoch, even Call A8927 appears at the red wing of the fea¬ 
ture (at ~8200 A in Fig. 1), whose presence disturbs the 
fitting. No other significant line of Call appears in our syn¬ 
thetic spectra between 4100 - 8500 A. 

The best-fit is made using two distinct line-forming re¬ 
gions during most of the pre-maximum epochs, but with a 
few exceptions. The PVF component does not seem neces¬ 
sary to fit the Call NIR feature at -10 days (see Fig. 1). 
Since the Call triplet is a very wide feature, the Call PVF 
component strongly overlaps with its HVF after -4 days. 
However, SYN-I—I- is not able to resolve such a complex spec¬ 
tral feature. When the optical depth profiles of the strong 


PVF and HVF components of Call starts to overlap in ve¬ 
locity space, SYN-k-1- assumes that only the stronger com¬ 
ponent at a particular wavelength (instead of the sum of 
the two) takes part in the line-forming process. This mod¬ 
elling assumption is completely different from the method of 
Childress et al. (2014) and S15, where the authors co-added 
the best-fit pseudo-absorptions of the PVF and HVF com¬ 
ponents (see below in this section). 

With the lack of observed Call H&K lines, we cannot 
create additional constraints on the Call HVF line forming 
region, and its properties become highly uncertain when the 
overlapping with the PVF component takes place. After - 
3 days, the PVF component shows a decent fit alone, and 
there is no direct sign of a resolved HVF component. Thus, 
we eliminated the Call HVF from our SYN-I—I- models after 
this epoch. Note that the ignorance of Call HVF does not 
necessarily imply the physical disappearance of the HVF line 
forming layer (see Sec 3.3). Similar result was presented by 
Vallely et al. (2015) for 2014J, where Call HVF was found 
absent in the near- and post-maximum SYNOW models. 

At -8 and -6 days the Vmin of the PVF component is 
slightly higher than Vphot, but later they run together (see 
Fig. 3). By extrapolating the two earliest Call PVF veloci¬ 
ties to -10 days, it is conceivable that the Call PVF starts 
as a DF shortly after the explosion, but because of its higher 
decline rate, later it decreases down to the photosphere, like 
other detached ions. The minimum velocity of the Call HVF 
starts at 30,400 km s“^, but declines fast until 4 days before 
maximum light (see Fig. 3). At the last two epochs (-4 and 
-3 days), when Call HVF is still detectable, the Vmin does 
not decline further, which may support the idea that Call 
HVF has a velocity floor in SN 2010kg (S15). 

In Fig. 3, the minimum velocities for Call and Sill are 
compared to the results for SN 2010kg of S15. Those au- 
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Figure 4. The velocities of ions with photospheric or detached line forming regions. The photospheric velocity is presented with the 
thick solid grey line, while the Si II HVF velocity (as the assumed upper limit of the SN ejecta) is plotted with the thin grey line. 


thors fit separate Gaussian profiles to the PVF and HVF 
components. They defined the PVF and HVF velocities as 
the Doppler shifts of the absorption minima of each Gaus¬ 
sian component, which are found to be consistent with the 
velocities from fitting P Gygni line profiles given by SYN-I—1-. 
As expected, the Ga ll velocities from these two studies are in 
good agreement. The difference between the HVF-velocities 
at the epoch of -10 days is probably due to the fact that 
we applied only the HVF component of Call for fitting the 
entire feature at this phase, while S15 used both the HVF 
and PVF components for fitting the same spectrum. 

4.2 Si II lines 

As for Can, the singly-ionized silicon shows both PVF and 
HVF. Although Sin forms another line at ~4900 A, this 
region is strongly overlapped by Fell features. As a result, 
the only Si ll feature where the HVF velocity is measurable is 
the A6355. The absorption components of the separate HVF 
-I- PVF can be seen by eye, just like the transition between 
the two components at 4 days before maximum light, after 
which the PVF becomes the dominant feature. Unlike the 
Call NIR triplet, there are no multiple lines of Sill in the 
region of A6355, thus, the HVF and PVF components of this 
feature stay resolved for longer time (Fig. 2). 

The velocity of the photospheric component decreases 
similarly to Vphot, but it is always higher by -^500 - 1500 
km s“^ than Vphot (see Fig. 3). These results show good 
agreement with the velocity values obtained by S15. The 
Vmin of Si II HVF does not change strongly, unlike the Ca ll 
HVF, the decline rate is only 200 km s“^ day“^ on the av¬ 
erage. After the B-band maximum, the Sill HVF does not 
disappear, contrary to most other SNe studied in recent lit¬ 
erature. Three days after maximum light its contribution to 


the absorption is still significant (see Fig. A2). At this time, 
its minimum velocity is ~18,000 km s“^, which is similar to 
the velocity floor of Call HVF in S15. Based on the high 
velocity floor of the HVFs in SN 2010kg, it seems possible 
that the origin of the HVFs might be related to some kind 
of interaction between the ejecta and an outer density en¬ 
hancement or GSM. 

4.3 Oxygen lines 

HVF of neutral oxygen is identified as the strongest O l fea¬ 
ture at A7773 (Parrent et al. 2012; Nugent et al. 2011), thus 
it seems reasonable to set both a PVF and a detached com¬ 
ponent in our models. At the early epochs, the contribution 
of OI PVF is negligible, thus we use only the DF compo¬ 
nent. This component shows the highest velocity among all 
the DFs in SN 2010kg (see Fig. 4), which is even higher 
than the HVF of Sill at -10 days. Later, the line forming 
region of O l sinks below the HVF line forming regions, but 
stays above any other ions, being at -^5000 km s~^ above 
the photosphere. Because of the continuous decreasing of its 
Vmin, we rather identify Ol DF than HVF. The behaviour 
of Ol suggests that the transition between HVFs and DFs 
may be continuous, and their line forming regions may not 
be separated sharply. 

At three days before maximum light, the O l DF can¬ 
not fit the red side of the A7773 feature alone, because the 
flux from the pseudo-emission of O l DF is too high to fit 
this spectral region. The application of Ol PVF improves 
the fitting till the end of the observed epochs. The veloc¬ 
ity of the PVF is bound to Vphot in our SYN-I—I- models 
(see Fig. 4). However, the two components (DF and PVF) 
cannot be seen separate, unlike in the case of Sill A6355. 
Thus, the detection of O l PVF is ambiguous in SN 2010kg. 
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Further testing of the simultaneous appearance of both the 
DF and the PVF components for O l may be possible in the 
NIR domain, where more oxygen features can be identified 
(Marion et al. 2009). 

Note that we also apply ionized oxygen (O ll) to achieve 
better fitting for the absorption feature at -^4400 A. We 
discuss the possible presence of O ll in the Section 5.3. 

4.4 Other ions 

Singly ionized magnesium (Mgll) is usually present in the 
early spectra of Type la SNe, either as PVF (Marion et al. 
2014) or DF (Parrent et al. 2012). Mgll may show more 
overlapping lines in the studied spectral range, but the red 
wing of the A4481 feature is relatively free of other blends, 
offering an excellent opportunity to measure its Vmin- Mgll 
is a typical detached ion in the atmosphere of SN 2010kg 
at the earliest epochs, but its initially high minimum veloc¬ 
ity decreases rapidly until it becomes permanently a PVF. 
This velocity profile is similar to the results from the anal¬ 
ysis of the SN 2014J (categorized as a broad-line la, just 
like SN 2010kg) optical spectra by Vallely et al. (2015) and 
Marion et al. (2015), where Vmin of Mgll declines through 
the epochs between -10 and +2 days. However, studying the 
near-infrared spectra of SN 2011fe (categorized as a core nor¬ 
mal subtype; Parrent et al. 2012), Hsiao et al. (2013) found 
that the Mg ll velocity has a floor at ~11,000 km s“^ between 
-10 and -1-10 days, which is supposedly due to a bottom limit 
of Mgll abundance enhancement in SN 2011fe. Since Mgll 
is produced only in carbon burning, the existence of Mgll 
Vmin floor suggests the presence of an inner edge of carbon 
burning in SN 201 Ife (Wheeler et al. 1998). The lack of the 
Mgll velocity floor in SN 2010kg implies that either carbon 
burning may reach deeper layers in broad line SNe la than 
in core normal SNe la, or, alternatively, mixing might be 
more effective in the former type of ejecta. 

Fe II and Fe ill usually appear together in the supernova 
ejecta, because they produce relatively high optical depth 
around T ~10,000 K (Hatano et al. 1999a). Both Fell and 
Fe III have several lines between 4100 and 5500 A, where the 
spectrum is strongly affected by Sill and Sll. Although the 
minimum velocity of the iron ions are ambiguous, it is com¬ 
mon to use Fell as a detached and Felll as a photospheric 
ion. In accord with these assumptions, we identify the lines 
of Felll as PVFs and Fell as DFs. The velocity of the latter 
one is ~20,000 km s“^ at the early epochs, which decreases 
linearly in time (see Fig. 4). Unlike Mgll, the Fell minimum 
velocity does not converge toward Vphot during the observed 
temporal evolution of SN 2010kg. 

The presence of Sll affects a wide range between 4200 
and 5200 A, but these lines are mostly weak. The two dom¬ 
inant absorption features are close to each other at A5463 
and A5641 A, often called the S ll ”W” feature. The bluer ab¬ 
sorption peak is overlapped by the emission part of the iron 
lines, but the redder peak is nearly free from other blending 
features. Because the line-forming region of S ll begins right 
above the photosphere (at least before maximum light), thus 
S II A5500 is ideal to test the Vphot values. 

To fit the absorption feature around ~5700 A, we 
adopted two more ions, Nal and Till, because Sill A5973 is 
not strong enough to fit the entire feature. The existence of 
both of these ions are possible in a Type la ejecta. The origin 


of titanium ions is usually explained by the thin helium-shell 
burning (Dessart & Hillier 2015; Valenti et al. 2013) in the 
literature, which is one of the possible triggering mechanisms 
of the Type la explosions. Ti ll is necessary only in the ear¬ 
liest spectra of SN 2010kg with minimum velocity values of 
20 - 22,000 km s“^. Our fitting uses Nal as a detached ion 
that shows nearly constant velocity (~13,000 km s“^) until 
maximum light. Note that the fit to the absorption feature 
around -^5800 A is very degenerate, thus, we cannot fully 
confirm the presence of either Nal or Till. 

As noted in Section 3.2, both Nill and Coll were also 
included in the initial models to reproduce the strong flux 
decrease below 4000 A as seen in every observed Type la 
SN spectrum (e.g. Foley et al. 2012a,b). This spectral re¬ 
gion was not covered by the HET spectra, which prevents 
constraining the optical depths of the Nill and Coll fea¬ 
tures. Therefore, we used these ions only to test whether 
their strong presence in the near-UV regime would affect the 
optical spectra. For this purpose we set the optical depths 
of Nill and Coll by hand, and varied them until the shape 
of the spectrum below 4000 A looked like those of other 
SNe la. The minimum velocities for these ions were fixed at 
Vphot- We found that the contribution of these ions is mi¬ 
nor, almost negligible in the optical spectrum above 4300 
A. Between 4100 and 4300 A a pseudo-emission peak from 
the Co II features may appear, but since our flux calibration 
is rather uncertain in this region, we cannot use this part 
of the spectrum to get realistic constraints on the strength 
of these features. Thus, we decided to omit both Nill and 
Coll from the final fits. Note that these spectrum models 
are valid only in the studied spectral range (4300 - 8500 A), 
and should not be treated as a complete description of the 
chemical composition of the ejecta of SN 2010kg. 

Although our SYN-I—I- models show a decent fit to the 
spectra of SN 2010kg, an unfit feature appears at -1 day in 
the observed spectra at ~5900 A, and it strengthens with 
time (see Fig. A2). Since there is no strong P Cygni feature 
in this area, whose pseudo-emission could explain the ob¬ 
served flux, it may originate from an emission source, which 
cannot be treated with SYN-I—I-. Such emission line is not ex¬ 
pected in Type la supernova around maximum light. Based 
on their theoretical model, Dessart et al. (2014a) showed 
that the forbidden transition of [Com] A5888 may be able 
to form such an emission line in delayed detonation Type la 
SNe. In their models the emission feature of [Com] A5888 
appears at a few days after maximum light, a week later 
than our best-fit SYN-I—I- models start to deviate from the 
observations at ~5900 A in SN 2010kg. 


5 DISCUSSION 

5.1 Temperature at the photosphere 

The Tphot values obtained from our SYN-I—I- models (~13 
- 15,000 K depending on epoch) are higher than the 
photospheric temperatures from some recent theoretical 
SN la models (e.g. Kasen et al. 2006; Jack et al. 2011; 
Dessart et al. 2014c; Kromer & Sim 2009), which are typi¬ 
cally under 11000 K. However, the lower temperature {Tphot 
< 11,000 K at maximum light) indicated by these models 
cannot fit the observed spectra of SN 2010kg, because the 
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Figure 5. The measured spectrum 5 days before maximum light (gray), the best fit synthetic spectrum with Tp^^t = 13,600 K (red), 
and a model spectrum having = 11,000 K (blue). The lower temperature model spectrum produces an inferior fit around 4500 

A. Above 6000 A the two models are indistinguishable. Note that the observed spectrum ends at 4100 A and the flux under 4300 A is 
uncertain. 


synthetic spectra do not reproduce the flux under 4700 A 
(see Fig. 5). This region may also be influenced by the un¬ 
known in-host reddening (see Sec. 2). Since we corrected the 
observed spectra only for the Milky Way reddening, these 
corrected flux values may be only lower limits of the true 
fluxes. Thus, the disagreement between the observed and 
the low-temperature model spectra could be even more pro¬ 
nounced than it is suggested by Fig. 5. 

Furthermore, the photospheric temperatures predicted 
by various, recently published SN la models, spread over a 
wide range, and some of them predict higher values. For ex¬ 
ample, Hachinger et al. (2013) suggested Tphot ~ 11 - 13,000 
K (depending on epoch) from the analysis of the SN 2010kg- 
like broad-line SN 2010jn. Dessart et al. (2014c) derived sim¬ 
ilar (T 12,000 K) temperature around maximum light in 
the velocity range of 10 - 16,000 km s“^, where most of the 
emergent spectrum is formed. Branch et al. (2005) applied 
Tphot ~ 13 - 14,000 K to fit the spectra of SN 1994D with 
SYNOW around and after maximum light. 

On the other hand, there are known cases when such 
temperatures do not fit the observed SN spectra. For exam¬ 
ple, Tanaka et al. (2008) modelled pre-maximum la spec¬ 
tra with a Monte-Carlo radiative transfer code, and found 
that Tphot ~ 17 - 18,000 K was necessary to get decent fits 
to spectra at 7 - 10 days pre-maximum. These are signifi¬ 
cantly higher temperatures than the ’’canonical” 11 - 13,000 
K ones mentioned above, even though Tanaka et al. (2008) 
found that such high temperatures are characteristics of 
low-velocity gradient (LVG) SNe, while high-velocity gra¬ 
dient (HVG) SNe, like SN 2010kg, tend to have lower Tphot- 
An extreme example for a SN la with high Tphot is SN 


2012dn (Chakradhari et al. 2014), where Tphot 20,000 K 
was found for the pre-maximum spectra. 

Another fact that further complicates the interpreta¬ 
tion of our derived Tphot values is that this parameter is 
very poorly constrained in a SYNH—h model, especially for 
Type la SNe where there is no well-defined continuum in 
the optical. Thus, Tphot given by SYN-I--I- cannot be taken 
at face value as a representative of a physically realistic tem¬ 
perature at the photosphere. Instead, it should be considered 
only as a fitting parameter, which describes the shape of the 
spectrum. 

Overall, it is concluded that even though our fits sug¬ 
gest higher Tphot values (~13 - 15,000 K) than what is com¬ 
monly adopted for SNe la (~11 - 13,000 K), these may not 
be directly related to the true temperatures around the pho¬ 
tosphere. Based on more realistic SNe la models (see above), 
such high temperatures may not be unexpected, though. 

5.2 Temperature distribution in the atmosphere 
of SN 2010kg 

Based on the calculations by Hatano et al. (1999a), we are 
able to create constraints on the excitation/ionization tem¬ 
perature distribution in the ejecta of SN 2010kg, considering 
on the optical depth values of the detected ions. As a first 
approximation, if we assume T^xc Tocai, then the pres¬ 
ence of photospheric Sill and the absence of Si ill suggests 
an upper limit for the temperature near the photosphere, 
because the turning point between their optical depth func¬ 
tions, where Si ill starts to dominate over Sill, is at ~11,000 
K. The optical depth functions of the sulphur ions also point 
toward T^xc Tocai 11,000 K. Similarly, the detection of 
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Rest wavelength (A) 


Figure 6. The absorption feature centered at ~4400 A in the -5 
days spectrum of SN 2010kg (grey line). The best-fit synthetic 
spectrum (including a C III detached line-forming region) is plot¬ 
ted with black solid line, the blue dashed line shows the synthe¬ 
sized spectrum without the C III contribution, while the orange 
dashed line shows the effect of replacing Cm with photospheric 
Si III. As seen. Si III produces an inferior fit compared to Cm. 



Rest wavelength (A) 

Figure 7. The measured spectrum 5 days before maximum light 
(gray), the best fit synthetic spectrum with Vp^^t = 14,200 km 
s“^ (black), and the effect of Vpf^^t = 11,000 km s“^ needed for 
the well-fitting Si III line (orange). Although the lower Vp^^t value 
improves the fitting of the 4400 A feature with Si III, the fit to 
the region around 5400 A gets worse, as also shown In the inset. 

doubly-ionized iron at photospheric velocity suggests a lower 
limit for Tiocai as ~12,000 K. 

The estimated excitation temperature regime is in con¬ 
flict with the Tphot values in our best-fit SYNH—h models. 
The reason for this discrepancy may be the same issue with 
the Tphot parameter that we discussed in Section 5.1. An¬ 
other probable reason is the lack of LTE in the atmosphere 
of SN 2010kg, as commonly found in other Type la SNe as 
well, which may make the excitation of different ions and 
species independent of the local temperature. Indeed, in a 
scattering-dominated SN ejecta, the color temperature of 
the emergent flux spectrum mimics the temperature of the 
deeper, thermalization layer, where true absorption starts to 
dominate over scattering (Jeffery & Branch 1990). Note that 
complete thermalization probably never occurs in Type la 


ejecta, thus, the emergent flux spectrum is not really Planck- 
ian (see e.g. Branch et al. 2005). As a result, T^xc = Tphot 
cannot be expected. 

Proceeding toward the outer regions in the atmosphere 
of SN 2010kg, the local temperature may be estimable based 
on the model T^xc values of the detached ions. The appear¬ 
ance of Fell at ~ 15 - 20,000 km s“^ (depending on epoch) 
and the absence of detached Felll suggests Toxc < 10,000 K 
(Hatano et al. 1999a). Based on our best-fit SYN-I—I- mod¬ 
els, Tiocai 8 - 9,000 K is most likely for this region of the 
ejecta, because the optical depth of Fe ll approaches 1 in this 
temperature regime. This estimate is also supported by the 
optical depth values of Ol, which stays as a detached ion 
through the pre-maximum epochs at velocities higher than 
Fell. 

The decrease of the local temperature from T 12,000 
K to 8 - 9,000 K between v 10,000 km s“^ and 20,000 km 
s~^ is in accord with the predictions from recent delayed- 
detonation models (e.g. Dessart et al. 2014c). 

5.3 The feature at 4400 A 

The identification of the narrow absorption-like feature at 
~4400 A rest wavelength turned out to be an interesting is¬ 
sue that is worth investigating. Many authors fit this feature 
with Sim A4560. Testing this assumption with the spectra 
of SN 2010kg, we find that Si ill is unable to fit the red wing 
of this feature, if its velocity is set equal to Vphot given by 
the fitting of Sll and Felll fines (see Fig. 6). In order to 
get a decent fit to this feature by Si ill the only possibil¬ 
ity is the decrease of Vphot by ^3000 km s“^. However, this 
has an impact on the whole spectrum and disturbs the fit¬ 
ting of the Felll A5156, the Sii W feature at ~5300 A and 
the pseudo-emission of Sill A6355 (see Fig. 7). We suggest 
that the absorption feature centred at ~4400 A is not likely 
formed by Sim, thus the existence of Sim in the ejecta is 
ambiguous. The only other feature which may support the 
presence of this ion (the small absorption notch at ~5500 A) 
is not strong enough to characterize its fine-forming region. 
It is still possible that the line at ~4400 A is caused by the 
mixed contribution of Sim and other ions, but because of 
its high uncertainty, we omit Si m from the fitting. 

On the other hand, it is found that this feature can 
be fit much better by using either C m, or O ll, as a de¬ 
tached ion (Fig. 8). Previously, the presence of Cm was 
identified in Type la spectrum only in the case of SN 1999aa 
(Garavini et al. 2004), although, they found Cm A4647 as 
a PVF. Chornock et al. (2006) also used C m to fit this fea¬ 
ture in the spectrum of the Type lax SN 2005hk, but they 
noted that this identification is rather ambiguous. 

For 2010kg, the minimum velocity for either O ll or C m 
must be '^4 - 5000 km s“^ higher than Vphot to get a decent 
fit. If real, this would put these ions among the fastest ones 
(see Fig. 4), but their line forming region stays definitely 
below that of the HVFs (see Fig. 3). This is the region where 
remnants of the unburned matter of the original C/O white 
dwarf are expected. 

Since either C m A4647 or O ii A4640 provides almost 
the same quality-of-fitting to the feature at 4400 A, we need 
to seek additional constraints to be able to identify this fea¬ 
ture. We calculate the optical depths for different ions of 
carbon (Fig. 9) and oxygen (Fig. 10) as functions of the lo- 
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Figure 8. The absorption feature centered at ~4400 A in the 
spectrum (grey line) of SN 2010kg, 5 days before B-max. The 
best fit synthetic spectrum (including a O II detached line-forming 
region) is plotted with the red solid line, while the black dashed 
line shows the synthesized spectrum with the C III contribution. 

cal temperature, similar to Hatano et al. (1999a). The grey 
shaded area represents the temperature range expected at 
the outermost ejecta region (Sec. 5.2). It is seen that Cm 
is unlikely to be really present in SN 2010kg, because this 
ion would require a local temperature well above 15,000 K. 
Instead, for temperatures lower than 10,000 K, it would be 
more feasible to find either C l or C ll, but there is no sign for 
any of these carbon ions in the optical spectra of SN 2010kg. 
All these pieces of evidence point toward the conclusion that 
the real presence of C ill is not probable. 

Although the ’’conventional” delayed-detonation mod¬ 
els (e.g. Kasen et al. 2006; Jack et al. 2011; Dessart et al. 
2014c) usually predict a decreasing temperature profile to¬ 
ward the outer ejecta, at least for the pre-maximum epochs, 
we also looked for alternatives, whether those might be able 
to support the appearance of C ill. For example, some of 
the pulsating delayed-detonation models by Dessart et al. 
(2014b) predict rising temperature profiles toward the outer 
region. Although the actual temperature values may depend 
on model assumptions, the range of the predicted outer tem¬ 
peratures is still not high enough to make the appearance of 
C III feasible. 

One other possibility might be the significant decrease 
of the electron number density with respect to the value of 

5 X 10® cm“^ adopted by Hatano et al. (1999a), in the outer 
ejecta. Solving the Boltzmann- and Saha-equations one may 
find that this would move the peak of the C ill ion optical 
depth toward lower temperatures. However, many orders of 
magnitude lower electron density would be necessary to get 
detectable amount of Cm around T ~ 10,000 K. Thus, we 
conclude that the presence of C ill cannot be justified. 

If, instead, we assume that the 4400 A feature is due to 
O II, then we find that the minimum velocity for O ll must 
be higher than Vphot by 4 - 5,000 km s“^ for every observed 
epoch. This is consistent with the O l velocities that are also 

6 - 7,000 km s“^ above the photospheric velocity (see Fig. 4). 
Since these two oxygen ions show the highest velocity from 
the ions forming DFs, they occur in the outer region of the 
SN 2010kg ejecta. Theoretical models (e.g. Seitenzahl et al. 
2013; Moll et al. 2014), which compute nuclear yield in Type 



Figure 9. The solid lines show the optical depth values of carbon 
ions with electron density Ae = 5 X 10® cm~^ and number density 
of carbon Ni = 5 X 10® cm“®, after Hatano et al. (1999a). The 
grey strip shows the local temperature range of the theoretical 
models in the outer part of Type la SN ejecta. 



Temperature (K) 


Figure 10. Optical depth values of oxygen ions with electron 
density Ne = 5 X 10® cm~^ and number density of oxygen Ni = 
5 X 10® cm“®. The grey strip shows the local temperature range 
of the theoretical models in the outer part of Type la SN ejecta. 


la supernovae, also show high oxygen abundance in the outer 
region of the SN la atmosphere. In these models, the oxy¬ 
gen appears with nearly constant mass fraction in the ejecta 
down to ~11 - 12,000 km s“^. 

In Fig. 10 we plot the optical depths of oxygen ions as 
functions of the temperature, again following the calcula¬ 
tions by Hatano et al. (1999a). It is seen that the predicted 
optical depth values for On, being between 0.01 and 0.1 in 
the grey-shaded temperature regime, are in good agreement 
with the values from our best-fit SYN-I--I- models (see also 
the tables in the Appendix). Moreover, this is also true for 
the OI optical depths found by SYNH—h- Thus, the simulta¬ 
neous presence of O l and O ll in the outer part of the ejecta 
of SN 2010kg seems to be plausible. 

Note that OI PVF is also detected in the spectra of 
SN 2010kg at 3 days before maximum light and afterward. 
This may raise the question: why does not the PVF com¬ 
ponent appear in the O ll feature as well at these epochs? 
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One possible reason might be the spatial separation of dif¬ 
ferent oxygen-rich line forming regions in the ejecta. Oxygen 
can occur in Type la atmosphere either as unburnt material 
of the WD or as the result of fusion process (Maeda et al. 
2010). Because of their spatial separation, oxygen PVFs are 
more likely formed in the freshly synthesized oxygen layer, 
which is expected to be located deeper in the ejecta. Thus, 
the number density of oxygen could be signihcantly differ¬ 
ent at different spatial regimes in the ejecta. If we assume 
lower amount of oxygen around the photosphere (< Vphot + 
2000 km s“^), only Ol may have high enough optical depth 
to produce a detectable optical line. This suggestion may 
be supported by the models of Seitenzahl et al. (2013) and 
Moll et al. (2014), where the computed oxygen abundance 
drops one or two orders of magnitude (depending on the 
model) below ~11,000 km s“^, but remains significant even 
at ~5,000 km s~^. 


6 CONCLUSIONS 

We present eleven observed spectra of the Type la super¬ 
nova SN 2010kg obtained between -10 and -|-5 days with 
respect to the B-band maximum, and their direct spectral 
analysis made via fitting synthetic spectra calculated with 
SYN-I--I-. The high cadence of the spectroscopic data give 
us the opportunity to reveal the velocity evolution of the 
identihed ions. We distinguish the observed absorption fea¬ 
tures into three groups: high velocity features (HVFs), pho- 
tospheric velocity features (PVFs) and detached features 
(DFs). The line-forming regions of HVFs show minimum 
velocities higher than 20,000 km s“^ at the earliest epochs, 
and about three days before maximum light they settle at 
a constant velocity between 18,000 and 20,000 km s“^. The 
typical HVFs are formed by Call and Sill in SN 2010kg. 
PVFs have (nearly) the same Vmin as the photosphere. The 
Vphot is usually determined by fitting the S ll and Fe ill lines, 
but Sill, Call, Ol, and Mgll also form PVFs at certain ob¬ 
served epochs. Detached ions have Vmin > Vphot but they 
are formed in the SN atmosphere, thus the minimum veloc¬ 
ity of DFs stay below the vhvf- Singly-ionized magnesium 
and iron ions start as detached ions at -10 days, but Mgll 
converges to the photosphere in time, while the Fell line 
forming region stays always above the photosphere. O l DF 
is detected with the highest velocity among the detached 
features. It might be formed in the unburned remnant of 
the progenitor white dwarf. 

The deep, narrow feature at -^4400 A is fit by Cm 
A4647 or On A4650 rather than Si ill A4560. This attempt is 
slightly in tension with other spectral analysis results pub¬ 
lished recently, thus, we study the possibility of doubly- 
ionized carbon and singly-ionized oxygen ions in the SN 
ejecta. Calculating the temperature dependency of the opti¬ 
cal depths for these ions, we find that high temperature (T 
~ > 16,000 K) is necessary to explain the presence of Cm 
in the outer regions of the ejecta. Since recent theoretical 
models suggest that the temperature decreases toward the 
outer atmosphere in a Type la SNe, the presence of C m is 
highly unlikely. 

Instead, the absorption feature at ~4400 A could be 
more likely due to O ll, which is able to produce a detectable 
optical feature at 7 - 9,000 K. This temperature regime is in 


good agreement with the temperature profile of SN 2010kg 
estimated from the optical depth functions of the detected 
ions. To our knowledge, this could be the first identification 
of ionized oxygen in a Type la SN atmosphere. 
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APPENDIX A: OBSERVED SPECTRA OF THE 
SN 2010KG, AND THE BEST FIT SYNTHETIC 
SPECTRA 

The following tables contain the SYN++ model parameters, 
which were calculated assuming exponential density profile 
and a reference velocity of 13,000 km s“^. Tphot is the tem¬ 
perature of the photosphere, Vphot is the velocity of the pho¬ 
tosphere, aO is the flux scaling parameter, while r is the 
optical depth of the ion, aux is the e-folding parameter of 
its distribution, Vmin is the minimum velocity of the its line 
forming region and Texc the Boltzmann-excitation temper¬ 
ature. 

This paper has been typeset from a hi® prepared by 

the author. 
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Figure Al. The observed spectra (grey lines) between epochs -10 and -4 days with respect to B-max, and the best fit synthetic spectra 
(red lines) 
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Figure A2. The observed spectra (grey lines) between epochs -3 and +5 days with respect to B-max, and the best fit synthetic spectra 
(red lines) 
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Table Al. The SYNH—h parameters of the best-fit synthetic spectrum 10 days before maximum light. The Vmin values 
s“^, while the T^xc values are in 10^ K. aO = 0.77, Vp^ot = 16.000 km s_i and Tp^ot = 15.000 K. 

are in 

10® km 

- 

Ca iipYF 

CaiipYF 

Si llpYF 

SillHYF OlpYF 

OloF 

On 

Mg II 

Sii 

Fell 

Fe III 

Nal 

Till 

log(tau) 

- 

4.20 

1.00 

3.35 

0.70 

-0.60 

0.58 

0.07 

0.30 

-0.10 

0.07 

-0.40 

’^min 

- 

30.4 

17.4 

21.5 

22.0 

20.7 

18.6 

16.0 

20.0 

16.3 

16.0 

17.4 

aux 

- 

5.50 

3.80 

1.67 

3.50 

1.62 

5.50 

2.35 

9.40 

1.40 

5.00 

3.00 

Texc 

- 

14.6 

10.5 

10.5 

5.0 

16.0 

6.0 

11.5 

5.0 

5.0 

10.0 

20.4 

Table A2. The SYNH—h parameters of the best-fit synthetic spectrum 8 days before maximum 

light. The Vmin values 

are in 

10® km 

s while the Texc values 

i are in 10^ 

K. aO = 0.76 

s 'Vphot = 15.200 km 

s_i and Tphat = 

14.600 K. 





- 

Ca iipYF 

CaiipYF 

Si llpYF 

SillHYF OlpYF 

OloF 

On 

Mg II 

Sii 

Fell 

Fe III 

Nal 

Till 

log(tau) 

1.00 

1.94 

1.59 

2.95 

1.74 

-0.90 

0.13 

0.95 

0.15 

-0.82 

0.77 

-0.74 

’^min 

17.5 

26.2 

15.8 

21.1 

21.1 

20.1 

16.0 

15.2 

19.3 

15.6 

15.2 

16.7 

aux 

6.00 

8.28 

1.46 

1.64 

2.00 

1.74 

5.84 

1.90 

8.41 

2.00 

2.29 

2.46 

Texc 

14.4 

14.4 

8.0 

8.0 

6.6 

16.0 

5.0 

9.6 

5.0 

5.3 

6.6 

20.5 

Table A3. The SYNH—h parameters of the best-fit synthetic spectrum 6 days before maximum 

light. The Vmin values 

are in 

10® km 

s while the Texc values 

i are in 10^ 

K. aO = 0.71 

> 'Vphot = 14.700 km 

s_i and Tphat = 

14.500 K. 





- 

C&llpYF 

Caiipv^p 

Si II PYF 

SillHVF OlpYF 

Oidf 

On 

Mg II 

Sii 

Fell 

Fe III 

Nal 

Ti II 

log(tau) 

1.26 

1.85 

1.37 

2.52 

0.77 

-1.07 

0.06 

0.80 

0.05 

-0.46 

0.50 

- 

'^min 

16.0 

23.4 

15.2 

20.9 

20.1 

18.7 

14.7 

14.7 

17.5 

14.7 

14.3 

- 

aux 

5.00 

4.80 

1.70 

1.70 

3.06 

1.56 

5.80 

9.50 

5.00 

9.95 

11.80 

- 

Texc 

15.5 

15.5 

7.2 

7.2 

6.3 

15.9 

13.0 

9.5 

5.0 

10.0 

11.8 

- 

Table A4. The SYNH—h parameters of the best-fit synthetic spectrum 5 days before maximum 

light. The Vmin values 

are in 

10® km 

s while the Texc values 

i are in 10^ 

K. aO = 0.68 

'Vphot = 14.000 km 

s_i and Tphot = 

14.300 K. 





- 

CaiipvF 

CsiUHVF 

Si llpYF 

SillHVF OlpYF 

O Idf 

On 

Mg II 

Sii 

Fell 

Felll 

Nal 

Till 

log(tau) 

1.50 

2.25 

1.46 

2.02 

0.28 

-1.26 

-0.12 

0.70 

0.27 

-0.36 

0.30 

- 

'^min 

14.0 

21.8 

14.6 

20.4 

19.3 

18.0 

14.2 

14.0 

17.0 

14.0 

14.0 

- 

aux 

3.90 

3.32 

1.47 

1.95 

4.80 

1.54 

6.20 

1.65 

4.30 

2.00 

2.30 

- 

Texc 

12.0 

12.0 

7.9 

7.9 

18.3 

16.0 

14.9 

8.8 

5.0 

10.2 

14.9 

- 

Table A5. The SYNH—h parameters of the best-fit synthetic spectrum 4 days before maximum 

light. The Vmin values 

are in 

10® km 

s while the Texc values 

are in 10® K. aO = 0.71 

> 'Vphot = 13.600 km 

s_i and Tphat = 

14.100 K. 





- 

Ca iipYF 

Ca iipY F 

Si II pyf 

SillHYF 

O loF 

On 

Mg II 

Sii 

Fell 

Felll 

Nal 

Ti II 

log(tau) 

1.85 

1.73 

1.37 

2.08 

0.10 

-1.58 

-0.15 

0.53 

0.20 

-0.60 

-0.05 

- 

'^min 

13.6 

20.4 

14.3 

20.2 

18.7 

17.5 

13.6 

13.6 

16.4 

13.6 

13.6 

- 

aux 

5.00 

4.85 

1.45 

1.77 

6.40 

1.60 

5.80 

1.80 

4.75 

1.90 

2.81 

- 

Texc 

12.9 

12.9 

6.4 

6.4 

25.0 

16.2 

24.4 

8.0 

5.0 

7.7 

5.0 

- 

Table A6. The SYNH—h parameters of the best-fit synthetic spectrum 3 days before maximum 

light. The Vmin values 

are in 

10® km 

s while the Texc values 

. are in 10^ 

K. aO = 0.68 

^ Vphot = 13.100 km 

s_i and Tphot = 

13.500 K. 





- 

Ca iipYF 

Ca iifiv F 

Si II pyf 

SillHYF 

O loF 

On 

Mg II 

Sii 

Fell 

Felll 

Nal 

Ti II 

logftau) 

1.45 

1.95 

1.11 

2.28 -0.15 

0.08 

-1.67 

0.20 

0.49 

0.16 

-0.83 

-0.12 

- 

'^min 

13.1 

20.0 

13.9 

19.8 13.1 

18.6 

16.8 

13.1 

13.1 

16.1 

13.1 

13.6 

- 

aux 

4.20 

3.40 

1.49 

1.50 1.29 

5.13 

1.55 

2.34 

1.74 

4.60 

2.64 

2.81 

- 

Texc 

11.7 

11.7 

6.4 

6.4 25.0 

25.0 

16.3 

6.8 

10.9 

5.0 

6.6 

10.2 

- 
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Table A7. The SYNH—h parameters of the best-fit synthetic spectrum 1 day before maximum light. The Vmin values are in 10^ km s 
while the Texc values are in 10^ K. aO = 0.66, Vp^ot = 12.200 km s_i and Tp^ot = 13.700 K. 


- 

Ca iipv^ F 

Caiipyp SiiipYF 

SiiiirvF 

O IpvF 

OloF 

On 

Mg II 

Sii 

Fell 

Felll 

Nal Till 

log(tau) 

1.80 

0.66 

1.74 

0.08 

-0.10 

-1.36 

-0.20 

0.17 

0.26 

-0.16 

-0.56 

I^min 

12.2 

12.8 

18.7 

12.2 

17.5 

16.1 

12.2 

12.2 

115.3 

12.2 

13.3 

aux 

4.12 

1.88 

1.63 

0.94 

5.00 

1.63 

2.10 

1.85 

2.00 

4.20 

3.38 

Texc 

9.8 

7.0 

7.0 

25.0 

25.0 

14.8 

10.6 

9.2 

5.0 

9.7 

12.0 


Table A8. The SYNH — h parameters of the best-fit synthetic spectrum at maximum light. 
Texc values are in 10^ K. aO = 0.68, Vp^ot = 11.600 km s_i and Tp^ot = 13.600 K. 

The V. 

min values are in 

10^ km 

s ^ , while the 

- 

Ca Tipy F 

Ca Hhvf 

SillpYF 

SiiiirvF 

O IpvF 

OloF 

On 

Mg II 

Sii 

Fell 

Felll 

Nal Till 

log(tau) 

1.80 

- 

0.54 

1.50 

-0.10 

-0.06 

-2.13 

-0.48 

-0.02 

0.15 

-0.48 

-0.90 

I^min 

11.6 

- 

12.2 

18.0 

11.6 

16.8 

15.70 

11.6 

11.6 

15.2 

11.6 

13.0 

aux 

4.00 

- 

1.92 

1.65 

1.70 

6.60 

1.61 

2.50 

1.80 

4.10 

0.50 

4.50 

Texc 

10.5 

- 

9.5 

9.5 

25.0 

25.0 

16.5 

20.0 

9.5 

5.0 

9.9 

22.0 


Table A9. The SYNH—h parameters of the best-fit synthetic spectrum 1 day after maximum light. The 
while the Texc values are in 10^ K. ctO = 0.71, Vp^ot = 11.000 km s—i and Tp^ot = 13.500 K. 

'^min values are 

in 10^ km s ^, 

- 

Caiipvp 

C&IIhVF SillpYF 

SillHVF 

O IpvF 

O idf 

On 

Mg II 

Sii 

Fell 

Felll 

Nal Till 

log(tau) 

3.15 

0.34 

1.34 

0.09 

-0.05 

-2.20 

-0.20 

-0.17 

0.05 

-0.56 

-0.60 

'^min 

11.0 

11.3 

18.0 

11.0 

15.6 

14.9 

11.0 

11.0 

14.5 

11.0 

13.0 

aux 

3.15 

3.00 

1.62 

1.30 

6.60 

2.20 

2.00 

1.87 

4.80 

0.81 

2.10 

Texc 

10.7 

7.6 

7.6 

25.0 

25.0 

15.6 

9.4 

8.7 

5.0 

10.0 

25.0 


Table AlO. The SYNH— h parameters of the best-fit synthetic spectrum 3 days after maximum light. The Vmin 
s“^, while the Texc values are in 10^ K. aO = 0.70, Vp^ot = 10.000 km s_i and Tp^ot = 13.000 K. 

values are in 10^ km 

- 

Caiipv^p 

Ca iifiVF 

SillpYF 

Si iihvf 

OlpVF 

Oidf 

On 

Mg II 

Sii 

Fell 

Feiii Nai Tin 

log(tau) 

2.03 

- 

0.09 

1.25 

-0.40 

-0.02 

-2.00 

-0.30 

-0.35 

0.00 

-1.60 

'^min 

10.0 

- 

11.1 

18.0 

10.0 

15.0 

14.3 

10.0 

10.0 

13.0 

10.0 

aux 

2.00 

- 

3.91 

1.42 

1.00 

6.00 

2.31 

1.00 

2.13 

6.38 

0.50 

Texc 

12.7 

- 

5.8 

5.8 

25.0 

25.0 

14.7 

5.0 

9.5 

5.4 

9.5 


Table All. The SYNH— h parameters of the best-fit synthetic spectrum 5 days after i 
s~^, while the Texc values are in 10^ K. aO = 0.76, Vp^ot = 9.400 km s_i and Tp^ot = 

maximum 

13.000 K. 

light. The Vmin 

values are in 10^ km 

- 

Ca Upv” F 

Ca Hhvf 

Si II PYF 

Si iihvf 

OlpvF 

OlpF 

On 

Mg II 

Sll 

Fell 

Feiii Nai Tin 

loff(tau) 

2.00 

- 

0.12 

1.15 

-0.06 

-0.29 

-1.98 

-0.30 

-0.75 

0.00 

-1.50 

'^min 

9.4 

- 

10.8 

17.5 

9.4 

14.4 

14.0 

9.4 

9.4 

11.9 

9.4 

aux 

3.90 

- 

3.11 

1.50 

2.50 

9.50 

3.20 

1.50 

0.85 

6.16 

0.66 

Texc 

10.6 

- 

6.0 

6.0 

25.0 

25.0 

14.2 

12.0 

10.0 

8.1 

9.8 
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